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1. INTRODUCTION 

The geomagnetic control of the earth's atmosphere is well under- 
stood. In the F-region and the topside ionosphere, non-electrical forces 

« 

transport plasma along the magnetic field lines only. In consequence, 
the world-wide distribution of ionization is strongly dependent on the 
dip angle. For example, the equatorial anomaly is roughly symmetrical 
about the dipole equator rather than the geographic. The same appears to 
be the case in the Jovian ionosphere (Mahajan, 1981). 

The influence of the magnetic field of Jupiter on its ionization 
pattern is one of several outstanding topics which need to be studied. 

Tan (1986) investigated the formation of equatorial anomaly in the Jovian 
ionosphere under a centered dipole model. Tan (1988) further studied the 
effect of the tilt of the Jovian dipole. The results were in broad 
agreement with those of a diffusive equilibrium model (Tan and Wu, 1981). 
In this study, an off-centered dipole model is constructed and its 
effects on the ionization pattern are investigated. \ 

1.1. THE MAGNETIC FIELD OF JUPITER 

Jupiter has the strongest magnetic field of any planet in the solar 
system.. Its magnetic dipole, not unlike the earth's, is tilted approxi- 
mately 10° from the rotational axis. In addition, significant quadrupole 
and octupole moments are present (Smith, et al., 1976; Acuna and Ness, 
1976). 

The Jovian magnetic field can be approximated as that due to an 
offset tilted dipole (OTD). A summary of the various OTD models can be 
found in Acuna, et al . (1983). The dipole is tilted approximately 10° 



in the meridional plane of longitude 220° in System III (Fig. 1). The 
dipole is located at approximately 0.1 Rj in the equatorial plane at 
longitude 150°, Rj being the radius of Jupiter. 

1.2. OBJECTIVE 

A qualitative inspection would reveal the fact that the tilt effect 
of the dipole would be maximum in the meridional planes of longitudes 
220° and 40° whereas the offset effect would be most prominent in the 
meridional planes of longitudes 150° and 330° (Tan and Wu, 1981; Tan, 
1988). The latter effect is presently studied by a simple offset dipole 
model without tilt. The dipole is assumed to be situated in the meridional 
plane of longitude Xjjj = 150°. Specifically, we study the differences 
in the ionization patterns in the meridional planes of Xjjj = 150° 

(Western sector) and Xjjj = 330° (Eastern sector) (see Dessler, 1983 for 
convention). 

In the Jovian topside ionosphere, photoionization production rates 
and the loss rates are extremely slow. Hence significant diurnal variation 
is not expected and a diffusive equilibrium model is justified. In this 
model, it would be appropriate to take the electron density at the base 
as the observed peak electron density, which displays the existence of 
equatorial anomaly (Mahajan, 1981). 



2. ECCENTRIC DIPOLE MODEL WITHOUT TILT 


The Jovian magnetic field can be approximated as that of a dipole 
situated in the equatorial plane of Jupiter at a distance a = 0.1 Rj from 
the center (cf. Acuna, et al . , 1983). In the present model, the tilt angle 
of the dipole from the rotational axis is neglected. Figure 2 depicts the 
meridional plane of longitude Ajjj = 150° which contains the dipole. The 
following systems of coordinates are defined in that plane: 

(1) Planetocentric cartesian (x‘, z'), 

(2) Planetocentric polar ( r 1 , 6'), 

(3) Dipole castesian (x, z), and 

(4) Dipole polar (r, 6) coordinates. 

The transformation from (r, 0) to (r‘, 0') is effected in the following 
stages: (r, 0) -A-*- (x, y) -B+ (x 1 , z') -C-*- (r 1 , 0'), where the transfor- 
mations A, B and C are given by 


A: 

X 

= r sin0 , 

(2.1) 


z 

= r cos0 , 

(2.2) 

B: 

x' 

= x + a , 

(2.3) 


z' 

= z , 

(2.4) 

C: 

r* 

= (x' 2 + z ,2 ) i , 

(2.5) 


0' 

= tan” 1 |-r • 

(2.6) 


Substituting Eqs. (2.1) - (2.4) in (2.5) and (2.6), we get 
r' = (r 2 + 2ar sin0 + a 2 ) 2 , 


(2.7) 
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and 0' = tan -1 (tane + -p sec0) . (2.8) 

The inverse transformations from (r 1 , 0') to (r, 0) are obtained by 
interchanging the primed and unprimed quantities and reversing the sign 
of a. 

r = (r* 2 - 2ar' sine 1 + a 2 )* , (2.9) 

and 6 = tan"^(tan0‘ -|rsec0') . (2.10) 

All equations (2.1) through (2.10) apply in the meridional plane of 
longitude Xjjj = 330° upon inversion of the sign of a. 

2.1. DETERMINING THE COORDINATES OF THE FOOT OF THE FIELD LINE 

For a dipole field line of definite L-value, the coordinates of the 
foot of the field line in the northern hemisphere (r^, 0 : rl, 0') can be 
calculated given the radial distance r o = R j + h 0 = R (say), where h Q is 
the altitude of the base of the ionosphere. For a diffusive equilibrium 
model of the Jovian ionosphere, h Q is approximately the height of the 
maximum electron density (800 km), whereas Rj = 71,200 km, whence R = 
72,000 km. 

The foot of the field line is the point of intersection of the 
circle r ,; = R and the dipole field line, whose equation is 

r = L Rj sin 2 0 . (2.11) 

Applying Eq. (2.11) in (2.7), we get 

L 2 R 2 sin 4 0 rt + 2aL R, sinV + a 2 - R 2 = 

Jo Jo 


o. 


( 2 . 12 ) 



One can solve for y = sin0 o numerically by the Newton-Raphson scheme 
(cf. Stark, 1970) 


fty) 


*i+i = 

■ w -j / 

y i f T TyTJ 

5 

(2.13) 

where f(y. ) 

■ L W + 

2aLRj y. 3 + a 2 - R 2 , 

(2.14) 

and f'(y.) 

■ 4L V^i 3 

+ 6aLRj y 3 , 

(2.15) 

whence 0„ = 
0 

sin _1 y , 


(2.16) 

r o 

LR, sin 2 0 rt , 
J 0 


(2.17) 

and 0^ = 

tan”^ (tan0 Q + 

7" sec V * 

(2.18) 


o 


The foot of the field line in the southern hemisphere then has the 
coordinates (r Q , 0 Q + 90°; r^, 0^ + 90°). 

2.2. DISTRIBUTION OF POINTS ALONG THE FIELD LINE 

A couple of transformations introduced by Kendall (1962) and by 
Sterling, et al . (1969) have been found convenient in controlling the 
distribution of points along the field line: 

2 4 

r Q cos© sin 0 cos0 

2 = 74 » 

r sin 0 

sinh(rg) 
sinh(rq 0 ) * 


and x 


( 2 . 20 ) 



Here is the value of q at the foot of the field line in the northern 
o 

hemisphere and r is an adjustable parameter which controls the distribution 
of points on the field line. A larger r for example distributes more 
points near the ends of the field line. 

The arc length s is related to the dimensionless quantity x through 


i_ = Y i_ 

3s " T 3x ’ 

(2.21) 

where y = - < 3 ^ V?* r cosh (rq) . 

(2.22) 


L r sin 0 sinh(rq Q ) 


2.3. LOCATING FIELD POINTS OVER THE LOCAL VERTICAL 

In order to study the vertical profile of the ionosphere over a 
locality at a certain latitude, it is advantageous to choose points on 
field lines over the local vertical. The coordinates of the field point 
(r 1 , 0'; r, e) are found from the intersection of the dipole field line 
with the line 0' = const. From Eqs. (2.10) and (2.11), we get 

2 3/2 

tan0 + -rf- - y -- tan0' = 0 . (2.23) 

LK J u 

Letting tan0 = u, Eq. (2.23) can be written as 

a n . ,,2x3/2 

f(u) = u + tt~ '■ 9 ] tan0' = 0 . (2.24) 

LR J u^ 


Then f'(u) 



3(1 + u 2 ) 3 / 2 
u 


2(1 + 


« bill 


u 


(2.25) 


and the Newton-Raphson scheme for u = tan© is given by 



Thus, given 0‘, the coordinates of the field point can be calculated. 


0 = tan _1 u , (2.27) 

r = LRj sin 2 0 , (2.28) 

and r' = (r 2 + 2ar sin0 + a 2 )* . (2.29) 

2.4. COMPONENT OF GRAVITY ALONG FIELD LINE 

The component of the gravitational acceleration g along the field 
line is given by (vide Fig. 3) 

g $ = - g si nl ' , (2.30) 

where I' is the true dip angle. If I is the dip angle in the centered 
dipole model, then 

I' = I + 0' - 0 . (2.31) 

Now, tanl = 2 cot© . (2.32) 

Substituting from Eqs. (2.32) and (2.8) in (2.31) and simplifying, we get 

, 2(tan0 + cot0) - % sec© 

I' = tan" 1 2 - . (2.33) 

sec^0 - -p sec0 (tan0 - 2 cot©) 

Taking into account the radial dependence of g, we have 

R 2 

g s = g 0 ^ sinl' , (2.34) 

where g 0 is the acceleration of gravity at the foot of the field line. 
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and r' and I' are given by Eqs. (2.7) and (2.33) respectively. 

In terms of a sole coordinate 9, 

r' 2 = L 2 Rj 2 sin 4 0 + 2aLRj sin 3 0 + a 2 , (2.35) 


and I* = tan 


2 (tan0 + cot0) - -pi— sec0 cosec 2 0 
-1 LR J 


2 

sec 0 


LR 


a 2 

sec0 cosec 0 (tan0 - 2 cot0) 


(2.36) 


2.5. COMPONENT OF CENTRIFUGAL ACCELERATION ALONG FIELD LINE 


The component of the centrifugal acceleration along the dipole 
field line is (vide Fig. 3) 


f = n 2 r' sine' cos (I'- e') . 


Since I'- 0' = I - 0, we have (from Tan, 1979) 


cos(I'- e 1 ) = 


3 cos9 sin9 
(3 cos 2 e + 1)' 


Further, from Fig. 2, we have 


r' sin0' = r sine + a , 


whence 


_ 3 fr(r sine + a) sine cose 
S (3 cos 2 e + 1 )* 


In terms of 6 only, Eq. (2.40) can be written as 


3 ft 2 (LRjSin 3 e + a) sine cose 
(3 cos 2 0 + 1 )* 


(2.37) 


(2.38) 


(2.39) 

(2.40) 


(2.41) 



3. THEORY 


The predominant ion in the Jovian ionosphere is proton (H + ). In 
the absence of electromagnetic drift, the governing equation of the ions 
along the magnetic field lines is given by (cf. Tan, 1986) 

If * fe( N V + BNV S fj(|) * o- BN . (3.1) 

where N = number density of ions/electrons, 

V s = velocity of ions along magnetic field lines, 

B = magnetic field intensity, 
s = arc length along B, 
t = local time, 
a = production rate of ions, 
and 8 = loss coefficient of ions. 

Now the topside ionospheric plasma can be approximated by a diffusive 
equilibrium distribution (cf. Bauer, 1969). This is even more true for 
the Jovian ionosphere, since the loss rate for protons is very slow. 
Under diffusive equilibrium conditions, Eq. (3.1) reduces to an ordinary 
differential equation 

35< nv s> + BNV s - 0 • < 3 ' 2 > 

In the Jovian ionosphere, the electrons, ions and neutrals are assumed 
to be in thermal equilibrium (Henry and McElroy, 1969; Nagy, et al., 
1976). The velocity V $ is then given by (cf. Tan, 1986) 

v s ■ 55 [rag s + " f s - 1 tff< NkT ” ' 


(3.3) 



where m = proton mass. 


v = ion-neutral collision frequency, 
k = Boltzmann's constant, 
and T = temperature. 

By substituting Eq.‘ (3.3) into Eq. (3.2), we have 


p , d_N + q. dN + R , N 
ds 2 ds 


= 0 , 


(3.4) 


with 


1 1 


2kT 
mv ’ 


(3.5) 


O' = 2k r 1 rfT + d_ f TN + I d_ f l n 
^ m *■ v ds ds v B ds 'B^ 


g +f 
~s s 

v 


(3.6) 


and R' = Hi r —(1 —1 + 1 -^1 R l-fln - ±J R d,JU 

m ^ ds^v ds' v ds ° ds'B'-* ds v ' v B ds'B' 


(3.7) 


Applying the transformation from s to x via Eq. (2.21), we get 


p^4 + q^+rn = o , 

dx^ dx 

(3.8) 

with P = y 2 P' , 

(3.9) 

9 * Q 1 . 

(3.10) 

and R = R‘ 

(3.11) 


Equation (3.8) can be rewritten as 



a dN 
a Sx 


+ bN = 


0 


> 


(3.12) 



where a = jj. , (3.13) 

and b = £ . (3.14) 

Equation (3.12) can now be solved by a standard numerical procedure. 
Rewriting N as y and dN/dx as z, we have two simultaneous first order 
differential equations 

^ = z = f^x, y, z) , (3.15) 

and if = " az “ ^ y» z ) • (3.16) 

Equations (3.15) and (3.16) can be solved by the Runge-Kutta scheme 
applied separately (cf. Stark, 1970) 

y i+1 = y n . + l (^ + 2k 2 + 2k 3 + k 4 ) , (3.17) 

z i+1 - z^ + -g (c^ + 2Cg 2Cg c^) , (3.18) 

where k-j = y-j* Zj)w » (3.19) 

C-| — ^2^ x i * 5 ^■j ^ s (3. 20) 

k c 

k 2 = f l (x i + I’ y i + T' Z i + T )vi * (3.21) 

k c 

c 2 = f 2 (x i + ?» y i + T~’ z i + r* )w * (3 - 22) 

k c 

k 3 = f l^ x i + ¥* y i + 2^’ z i + ’ (3.23) 

k c 

c 3 = V x i + 7» yi + jp z i + jr) w » (3.24) 



f-| (x^+ w, y^+ k 3 , z^+ c 3 )w , 


(3.25) 


and c 4 = f 2 (x 1 -+ w ’ y i + k 3* z i + c 3^ w 


(3.26) 


In Eqs. (3.19) through (3.26), w is the step size, i.e., Ax. 


3.1. STATIC DIFFUSIVE EQUILIBRIUM 

The condition of static diffusive equilibrium, which is applicable 
to the topside ionosphere, is given by V $ = 0 (cf. Bauer, 1969). Then Eq. 
(3.3) can be integrated, by separation of variables, to give 


N = 


NT 
o o 


4 IFr ( 9s + V ds 


The differential arc length in terms of 0 is 


(3.27) 


ds 


LRj (3 cos 2 9 + 1)* sine de 


(3.28) 


Equation (3.28) can be integrated from the foot of the field line in the 
northern hemisphere (cf. Chapman and Sugiura, 1956) to give 


s = LRj [cos0 Q (3 cos 2 0 q + 1)* - cos0 (3 cos 2 0 + 1)^] 



/3 cos6 Q + (3 cos 2 0 q + 1)^ 
/3 cose + (3 cos 2 0 + 1)* 


(3.29) 


In terms of the variable x (whose values range from -1 at the foot 
of the field line in the northern hemisphere through naught at the dipole 
equator to +1 at the foot of the field line in the southern hemisphere), 
Eq. (3.27) assumes the form 
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N = 


t x m , , - x dx 

NT / - 1 2kT t9s V T 
0 0 „ 1 
— t- e 


where y is given by Eq. (2.22) 


(3.30) 


3.2. INPUT DATA 

In the diffusive equilibrium model, it is appropriate to take the 
electron density at the base as the peak electron density. The observed 
peak electron density in the Jovian ionosphere, displaying equatorial 
anomaly, has been fitted with the expression (Tan, 1986) 

N Q = 0.18 + 0.0367 X 3/2 - 0.004056 X 2 , (3.31) 

where X is the absolute value of the latitude in degrees and N Q is 
5 3 

expressed in 10 electrons/cm . Figure 4 shows the latitudinal variation 
of N q as given by the expression (3.30). 

The altitude of the base is taken as the height of the observed 
peak electron density, which is approximately equal to 800 km (see, for 
example, Mahajan, 1987). 

The altitudinal variation of the temperature is taken from Atreya, 
et al. (1981) and Festou, et al . (1981). The temperature is assumed to 
increase linearly from 600°K at 800 km to 1200°K at 1400 km. The 
temperature above 1400 km is assumed to have the constant value of 1200°K. 
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4. RESULTS AND DISCUSSIONS 

Figures 5 through 9 illustrate the results of our calculations. In 
Fig. 5, contour plots of the electron density in the longitudinal planes 
of both Western and Eastern Sectors are shown together with the dipole 
field lines. In general, the electron densities are higher in the Western 
Sector than in the Eastern. This is a direct consequence of the con- 
figuration of the dipole field lines in the two Sectors. In the Western 
Sector, the ions can diffuse along the field lines to higher altitudes 
from the base at the same latitude, thereby producing greater electron 
densities in the pi asmasphere. 

Also shown in Fig. 5 are the loci of points at which the components 
of gravitational and centrifugal forces along the field lines are equal 
and opposite. In the region of centrifugal dominance, the electron density 
slowly increases with distance, thus creating a region of enhanced 
electron density away from the planet. As a result, the troughs of 
electron densities are found at midlatitude heights (A = 20° to 30° in 
the Western Sector and A = 20° to 25° in the Eastern Sector) rather than 
above the equator. Above the equator, the electron density profile 
possesses a minimum which is the saddle point in the electron density 
contour. 

Figures 6 and 7 show the vartical profiles of electron density in 
the two Sectors above the equator (latitude A= 0°) and a midlatitude 
location (A = 30°). As stated earlier, the latitudinal variation of the 
electron density at the base was given Eq. (3.31) to reflect the 
observed electron density. The electron densities were generally lower 
above the equator than at the midlatitude location owing to the boundary 



conditions used. However, this situation is reversed above 6500 km in 
the Western Sector and above 9500 km in the Eastern. This means that the 
equatorial anomaly does not translate into the plasmaspheric heights. As 
stated above, the electron densities are always greater in the Western 
Sector compared with those in the Eastern. 

The plasma scale heights were also greater above the equator and 
in the Western Sector. The difference in the scale heights between the 
two Sectors above the equator (1700 km vs. 1360 km) was particularly 
significant. This difference ought to betray itself if enough 
observational data were available. The calculated scale heights were 
generally larger than those observed which ranged from 540 km to 1040 km 
(cf. Mahajan, 1987). 

Figures 8 and 9 show the vertical profiles of electron density in 
the plasmasphere above the equator and the midlatitude location 
respectively. The gradients of the curves become steeper with altitude 
and consequently the scale height becomes larger. This is due to the fact 
that the centrifugal force on the corotating plasma increases with the 
radial distance as the gravitational force decreases. 

As the electron density curve becomes steeper with altitude, the 
electron density may undergo a minimum and actually increase with 
altitude thereafter. Figure 10, which is an extension of Fig. 8 to 
include higher altitudes actually shows this. This would be true so long 
as the assumption of the corotation of plasma is valid. As for the limits 
of corotation on Jupiter, the reader is referred to an article by Hill 
(1979). 



4.1. FURTHER STUDY 


This simplified model has been conveniently applied to study the 
ionization patterns in the meridional planes of maximum offset effect. 

The ionization patterns in the other meridional planes would be inter- 
mediate between the two extreme cases studied. However, in the other 
sectors, the mathematical expressions would be much more complicated. 
Further studies on the effect of the planetary magnetic field on the 
ionosphere would include the following topics. 

(1) Saturnian ionosphere model in the offset dipole approximation. The 
magnetic field of Saturn is more purely dipolar than Jupiter's with the 
dipole axis within 1° of the rotational axis. The dipole is off-centered 
by about 0.04 R$ to the north. Here the longitudinal differences in the 
ionization patterns would not be significant. But the offset effect would 
produce differences between the ionizations in the northern and southern 
hemispheres. A simplified model along the lines of this study can be 
easily constructed. 

(2) Offset dipole model with tilt. The most general dipole model is that 
of an offset and tilted dipole. For Jupiter, the dipole could be taken 
as lying in the equatorial plane. This assumption fetched a great deal 
simplification in the mathematical expressions. However, in the case of 
Uranus or Neptune, this luxury cannot be afforded. The magnetic dipoles 
of both these outer planets are highly eccentric and tilted. Here one 
needs a general offset and tilted dipole model to study the magnetic 
field effect on the ionization patterns. Needless to say, this will be 

a much formidable challenge for the modeller. 
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Fig. 2. Planetocentric Cartesian (x',z')» Planetocentric Polar (r 1 , 6'), 
Dipole Cartesian (x, z) and .Dipole Polar (r, 0) Coordinates. 
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The Jovian Plasmasphere in Offset Dipole Model* 

A Tan (Department of Physics, Alabama A&M University, 
Normal, AL 35762; 205-851-5313) 

The magnetic field of Jupiter can be approximated by 
that of a dipole situated in the equatorial plane at a 
distance of O.lRj from the center at longitude 150° in 
System III. The offset effect would be most prominent 
in the meridional planes of longitudes 150° (Western 
Sector) and 330° (Eastern Sector) . A diffusive 
equilibrium model in an offset dipole approximation is 
constructed to study the ionization patterns in the 
two sectors. The electron density is generally greater 
in the Western Sector than in the Eastern. The plasma 
.scale height is also higher in the Western Sector. The 
plasma scale height is greater at the equator than at 
a midlatitude location (A = 30°). Even if the peak 
electron density is assumed lower at the equator to 
reflect equatorial anomaly, at plasmaspheric heights, 
the electron density above the equator exceeds that at 
midlatitide regions, owing to the centrifugal 
dispersion of plasma there. 

*This study was supported by NASA Grant NAGW-1067. 
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